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Station (WES). The study was performed, and this report prepared, by Messrs. 
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Division. The WWD and its personnel were transferred to the Coastal Engineer- 
ing Research Center (CFRUC) of WES on 1 July 1983 under the supervision of 
Dr. Whalin, Chief of CERC. 

This report describes the application of numerical hydrodynamic models 
to Pamlico Sound and adjacent areas. Following model calibrations and 
verifications, ar. evaluation of proposed jetty configurations for Oregon 
Inlet was made. All numerical computations were performed on the CRAY 1 
and Control Data Cyber computers of Boeing Computer Services in Bellevue, 
Washington. 

Commanders and Directors of WES during this investigation and the rrep- 
aration and publication of this report were COL Nelson P. Conover, CE, and 


COL Tilford C. Creel, CE. Technical Director was Mr. F. R. Brown. 
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NUMERICAL SIMULATION OF OREGON I.LET CONTROL STRUCTURES’ 
EFFECTS ON STORM AND TIDE ELEVATIONS IN PAMLICO SOUND 


PART I: INTRODUCTION 


Background and Objectives 


1. North Carclina's Atlantic coast has two large lagoons, Pamlico 
Sound auu Albemarle Sound, separated from the ocean by a strand of barrier 
islands known as the Outer Banks (Figure 1). The OCuter Banks stretch approxi- 
mately 190* miles between Virginia Beach, Virginia; Cape Hatteras, and 
Merehead City, North Carolina. Only three openings in the barrier islands 
provide navigational access between Pamlico Sound and the Atlantic: Hatteras, 
Ocracoke, and Oregon Inlet. Oregon Inlet, the northernmost of the three 
openings, provides the fishermen and boaters of Roanoke Island and Albemarle 
Sound the only practical means of reaching the ocean. Since the creation of 
the present~day Oregon Inlet by a hurricane in 1846, navigation there has 
remained hazardous due to channel mcvements and shoals. Location of the 
channel] also nas meved steadily south du2 to sediment transport by alongshore 
currents in the past 100 years. 

2. In 1970, Congress authorized the Manteo (Shallowbag) Bay project 
which included the stabilization of Oregon Inlet with a dual jetty system. 
Since the addition of jetties will change the hydrodynamics of the inlet, 
proposed structures need to be evaluated for their impact on normal tidal 
conditions and storm curges. 

3. This study evaluates several proposed jetty configurations for 
Oregon Inlet with numerical hydrodynamic models of the entire lagoon system 
and the area around the inlet. The models are calibrated by computing tides 
and acjusting model parameters until the computations match observed tides. 
The calibrations are verified by storm surge simulations of two historical 
storms (Hurricane Donna in 1960 and the March 1962 northeaster) that affected 
the inlet. These calculations are showu to agree with marigrams recorded dur- 


ing the storms at several tide stations on the Atlantic coast and in the bays. 


* A table of factors for converting U. S. customary units of measurement to 


metric (SI) units is presented on page 3. 


After verification, the effects of the jetties are estimated with the models 


for both tide and storm conditions. 


Approach 


4, The numerical hydrodynamic code used in this study is WIFM, the WES 
Implicit Flooding Model (Butler, in prenaration). WIFM employs finite dif- 
ference methods to approximate the vertically integrated Navier-Stokes equa- 
tions. Several special features of WIFM are important to this study. The 
flooding and drying of tidal flats and low-lying lands are accurately simulated 
in WIFM. Variably spaced finite difference grids are used to maximize accuracy 
in the hydrodynamic simulations, without increasing computer costs. WIFM can 
approximate jetties, and other small topographic details, as thin barriers 
between computatienal cells. The implicit formulation used in WIFM's solution 
scheme allows for larger time-steps than do explicit solution schemes. 

5. The hydrodynamic features of the Carolina coast and Oregon Inlet 

are simulated to three levels of accuracy in this study by the use cf three 
different finite difference grids that cover different portions of the study 
area (Figure 2}. The three grids can be linked together, wherein an internal 
boundary in a larger grid can provide boundary conditions to the next small- 
est grid. The offshore grid (Figure 3) approximates the entire Carolina 
coast with 3,186 cells, making computaticns with it inexpensive. This grid 
extends seaward to the continental shelf in order to model the effects of 
long-wave shoaling caused by the shelf. Oregon Inlet is only approximated by 
one cell in the offshore grid, so computations from this model can only 
provide rough estimates of hydrodynamics in the sounds. The nearshore grid 
(Figure 4) models the Outer Banks in much finer detail with 9,009 cells. 
Tides and surges can be accurately simulated throughout the bays and inlets 
with this grid. Oregon Inlet is modeled with enough detail in this grid to 
allow for calculations with crude jetty approximations. This grid provides 
all boundary conditions to the most detailed of the three grids used in this 
study, the shore process grid. 

6. The shore process grid (Figure 5) covers less than 85 square miles 
around Oregon Inlet with 4,620 computational cells. This model provides high 
resolution at the inlet and in the surf zone. The shore process model predicts 


circulation patterns at the inlet, predicts effects of jetties on inlet flow, 
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Figure 5. Shore process model computational grid 
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and provides flow data for the sediment transport studies of Houston (in prep- 
aration). The e:fects of cuannel dredging can also be studied with this 
model. 

7. Scorcm simulations in WIFM require the input of wind stresses and 
barometric pressures into the computations. The wind and pressure fields for 
Hurricane Donna are simulaced with the Standard Project Hurricane (SPH) model 
(Graham and Nunn 1959; National Oceanic and Atmospheric Administration 
(NOAA) 1972; Schwerdt, Ho, and Watkins 1979). SPH is a parametric model that 
characterizes a hurricane by its radius to maximum winds, central pressure 
deficit, forward speed, and other variables. The parameters needed to re- 
Produce a historical storm can be cclilected from surface weather charts. The 
winds and pressures of the March 1962 ncrtheaster are estimated from digitized 
daca provided by the Wave Information Study (WIS) at the U. S. Army Engineer 
Waterways Experiment Station (WES) (Corson, Resio, and Vincent 1980; Resio, 
Vincent, ard Corson 1982). WIS data are referenced to an earth-coordinate 
system, and a computer code was developed to interpolate this information to 
WIFM grids. Wind velocities and pressures for each WIFM cell are read in 
during the hydrodynamic calculations, and WIFM converts wind velocities in 


surface stress by Charnock's relation (Garrett 1977). 
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PART II: COMPUTATIONAL TECHNIQUES 


Equations of Motion 


8. The hydrodynamics of the numerical model WIFM are derived from tne 


Navier-Stokes equations in a Cartesian coordinate system (Figure 6). The 


WATER SURFACE 


— BENCHMARK DATUM 


Figure 6. Coordinate system for problem formulation 


long wave approximations of small vertical accelerations and a homogenous 


fluid yield vertically integrated two-dimensional equations of continuity and 


momentum: 
Continuity 
CLAS: a) es 
TES (ud) + 3y (vd) = R (1) 
Momentum 
; 1/2 
Ju du du C gu Z 2 
— ais at a ; y—- = 
el ae ay eee oe I 
Cd 
2 2 
au dou 
- € 72 + 372 + Fy = 0 (2) 
x y 
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dv ov 3 3 2 
< << — = ego + 
svar Oar es Wilinee sh fu + g By (n n) 7) (u Vine) 
Cd 
ae av 
Se ae te eee (3) 
ae we NEOs y 


The dependent variables in the problem are 1 ,Uu , and v_, which represent 
the water-surface elevation above datum and the vertically averaged water 
velocities in the x- ond y-directions, respectively.* The other variables in 
the equations are: h , the local ground (cell) elevation above datum; 
d=n-h, the total water depth; t , time; f , the Coriolis parameter; 

C , tne Chezy coefficient ror bottow friction; ¢ , the eddy viscosity coet- 
ficient; g , the acceleration due to gravity; and R , the rate of water 
volume change in the system due to rainfall or evaporation. The coefficient 
me accounts for hydrostatic water elevations due to atmospheric pressure 
differences, and ES and a are terms representing external forces such as 
wind stress. 

9. The computational grid used for the finite difference approximations 
in WIFM employs a stretching transformation for each of the two coordinate 
directions (x and y). This transformation from physical space to computational 
space offers the major advantage of allowing increased grid resolution in 
areas of interest, by controlling the arbitrary constants a, b, and c in 


the equation: 


x =a + ba (4) 


Physical distances are defined by x , and the computational grid lines ir 
each direction are defined by positive integer values of a. The values of 
a,b, and c are determined not only from desired grid resolutions, but 
also from the requirement that the derivative a be continuous everywhere. 
Many stability problems commonly associsted with variable grid schemes are 


eliminated due to this continuity constraint. The transformation is applied 


—_—: a 


* For convenience, symbols and unusual abbreviations are listed and defined 
in the Notation (Appendix A). 
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aaa 
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to each coordinate direction independently in order to maximize grid resolu- 
tion in areas of hydrodynamic importance and to minimize computational cells 


in the far field. 


Numerical Method 


10. The alternating-direction-implicit (ADI) method has been used by 
Leendertse (1970) and others to solve the two-dimensional equations of motion. 
When the advective terms are included in the momentum equations (Equations 1-3) 
the ADI method has encountered stability problems. Weare (1976) indicates 
that the problem arises from approximating advective terms with one-sided 
differences in time, and suggests the use of a centered scheme with three 
time-levels. WIFM employs a centered stabilizing-correction (SC) scheme which 
is second-order accurate in space and time, and boundary conditions can be 
formulated to the same order accuracy. Details of the SC scheme can be found 
in Butler (in preparation) and a general development is presented in the 
following paragraphs. 

11. The linearized equations of motion can be written in matrix form 


as: 
U. + AU, + Bu = 0 (5) 


where 


The SC scheme for solving Equation 5 is 


k-1 
k = as = 
(1 + ee) U (1 hy ee) U (6) 
Gp anos = peat gest (7) 
yi y 
where 
1 At 1 At 
ry 5 Ae Aé and A 2 ay hss 
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The quantities oe and Oe are centered difference operators and the super- 
script k counts time-levels. The starred quantities can be considered 
intermediate values for variables at the (k+l) time-level. 

12. The first step in the SC procedure computationally sweeps the grid 
in the x-direction, with the second step sweeping the y-direction. Completing 
both sweeps constitutes a full time-step, advancing the solution from the a 


time-level to the (k+l) time-level. The form of the difference equations 


for the x-sweep is given by 


Fie IC eriaeete, Cutdta 1d) +5 6) (Say a2 ee? 
se ut 4) + cK nt t 4) = 0 (9) 
Tee (we - VI) + Bg nh?) = 0 (10) 

and the y-sweep by 
zl ttle sy oy ral eenbs ees he (12) 
re gk (12) 
oo pee n # : (Atl Kalyd ig (13) 


13. Noting that v* in Equation 10 is only a function of previously 
computed variables at the (k-1) time-level, its substitution into Equa- 
tion 13 and the substitution of u* (Equation 12) into Equations 8 and 9 
yield the simplified forms 


x~sweep 
1 k-1 1 k+1 -1 1 k-1 & 2 
a ee ae este, ate = 
Tat (n n )+ Phx 5 (u- “d+u “d) + ; oa (v ~d) ) (14) 
Mle hake) aut k=) g k-1, _ 
ThE (u -u §) + aan 8 (n® +n) = 0 (15) 


y-sweep 


1 k+1 1 k+1 k-1 
ancien — n* aS - = 0 (16) 
TAC (n n*) + Shy . (v d Vv d) 
le iH wit) ee es C “ ne 4) ZG (17) 
2At 2dy sy 


14. The details of applying the SC scheme to Equations 1-3 can be found 
in a report by Butler (in preparation). The diffusion terms of Equations 2 
and 3 are also represented with time-centered approximations. The inclusion 
of diffusion terms contributes to the numerical stability of the scheme 
(Vreugdenhill 1973), and serves to somewhat account for turbulent momentum 
dissipation at the larger scales. While the resulting finite difference forms 
of Equations 1-3 appear cumbersome, they are efficient tu solve. Application 
of the appropriate equation to one row or column of the grid (the "sweeping" 
process) results in a system of linear algebraic equations whose coefficient 
Matrix is tridiagonal. Tridiagonal matrix problems can be solved directly, 
without the cost and effort of matrix inversion. 

15. The computational time-step for the SC scheme in WIFM js largely 
governed by simpie mass and momentum conservation principles. The maximum 


time-step for a problem is characterized by: 


nAs 


At => 


(18) 
where V is the largest flow velocity to be encountered at a cell with its 
smallest side length AS . The parameter n is o: order 1. So, the time 
step is constrained by the smallest cell width which contains the highest flow 
velocity. In physical terms, Equation 18 requires that the flow cannot move 


substantially farther than one cell width in one time-step. 


Boundary Conditions 


16. WIFM allows a variety of boundary conditions to be specified, which 
can be classified into three groups: open boundaries, land-water boundaries, 
and thin-wall barriers. 


Open boundaries 
17. When the edge of the computational grid is defined as water, such 


16. 


TT ae ee - 
A SREY SOSA. ES 


a aT aso alia 


as a seaward boundary or a channel exiting the grid, either the water eleva- 
tion or the flow velocities can be specified as an cpen boundary-condition. 
This information can be input to WIFM as tabular data, or constituent tides 
can be calculated within the code during the time-stepping process. Open 
boundaries for a grid can be saved from a specitied internal boundary of 
another griu so that computational grids can be linked together. Grid linking 
is used in this study in order to model large coastal areas inexpensively and 
to supply correct boundary conditions to the shore-processes grid. 
Land-water boundaries 

18. WIFM allcws land-water boundaries to be either fixed or variable to 
account for flooding in low-lying terrain. Fixed boundaries specify a no-flow 
condition at the cell face between land and water. The position of a variable 
boundary is determined by the relationship of the water elevation at a "wet" 
cell to the land elevation at a neighboring "dry" cell. Once a water eleva- 
tion rises above the level of adjacent land height, water is initially moved 
onto the "dry" cell by using a broad-crested weir formula (Reid and Bodine 
1968). When the water level on the dry cell exceeds some small value, the 
boundary face is treated as open and computations for ) , u , and v are 
made at the now "wet" cell. Drying is the inverse process, and mass is con- 
served in these procedures. 
Thin-wall barriers 

19. These barriers are defined along cell faces and are of three 
types: exposed, submerged, and overtopping. Exposed barriers allow no flow 
across a cell face. Submerged barriers control flow across a cell face by 
using a time-dependent friction coefficient. Overtopping barriers are dynamic: 
they can be completely exposed, completely submerged, or they can act as 
broad-crested weirs. The barrier character is determined by its height and the 


water elevations in the two adjoining cells. 


Grid Connections 


20. Application of the embedded grid concept employed in ths study 
requires a transter of hydrodynamic information frem one grid to another. 
Specifically, data are transferred from the offshore grid to the nearshore 
grid and thence to the shore process grid. In practice, the shore process 


mode]. was always driven by a nearshore model simulation, and the offshore 
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model was used to drive the system only for storm events. 

21. The two embedded grids utilized in this study (nearshore and shore 
process) contain successively finer resolution with resultant greater detail. 
Therefore stability requirements as governed by Equation 18 necessitate the 
use of successively smaller time-steps in simulations with these models. Data 
are transferred between coupled grids at the outer boundary cells of the 
embedded grid. Coupling grids are designed such that the cell size Ax and 
the derivative as (in the direction orthogonal to the coupling boundary) are 
equivaient in the coupling cells of both the driving and driven grids. In the 
direction normal to the coupling boundary, cell expansion and contraction are 
allowed to proceed independently, thus requiring spatial interpolation of the 
transferred data. A temporal interpolation is also required due to the small 
time-step needed by the embedded grid. Care must be taken to ensure that 
coupling grids contain similar water volumes within the embedded grid cover- 
age, particularly at the coupling boundaries in order to preserve hydrodynamics 
between the grids. 

22. The data transferred between models consist of nn, u , and v at 
each coupling cell. A utility program performs the necessary spatial and 
temporal interpolations between coupling grids. It should be noted that 
connections are not totally dynamic in the sense that the simulations are 
independent ard not concurrent. Communication is unidirectional from the 


driving to the driven model. 
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PART II?: PROTOTYPE DATA 


23. Field meisurement data required for model calibrations and verifi- 
cation were taken from existing data bases. No new collection efforts were 
undertaken during the course of this study. The following is a summary of the 


prototype data used in the investigation and its sources. 


Tidal Data 


24. A hydrautic model study to determine effects of Oregon Inlet 
stabilization by jetties was conducted at WES and reported by Hollyfieid, 
McCoy, and Seavergh (1983). Representatives ot the National Ocean Survey 
(NOS), U. S. Army Engineer District, Wilmington (SAW), and WES tormulated a 
data collection program. Field surveys were carried out in the spring and 
summer of 1975. Efforts included collecting tide and current data, primarily 
in the vicinity of Oregon Inlet, ana prototype hydrography/topography data 
within the area reproduced in the hydraulic model. The tide data were sub- 
sequently analyzed for their tidal constituents. Constituent data form allows 
for isolating the astronomical event and was used both to force the WIFM model 
and for comparison with model results. Tide elevation data outside the Oregon 
Inlet area were obtained (in constituent form) from NOS. Figure 7 (pag2 24) 
shows locations of field stations in the offshcre grid region. Since the My 
constituent contains about 90 percent of the tidal energy in the Pamlico Sound 
area, calibration computations were simplified by considering only an M, tide. 


2 


Table 1 delineates the M, amplitude and epoch (relative to Greenwich) for each 


field station. 


Storm Surge Data 


25. Two historical storms were modeled in this study: Hurricane Donna 
in 1960 and the extratropical March 1962 storm. Meteorological data for 
Hurricane Donna were obtained from NOS. These jncluded a list of storm param- 
eters, track definition, and surface windfield analysis. Table 2 presents a 
time variation of these storm parameters. The coincident predicted astro- 
nomical tide was reconstructed from constituent tidal data (particularly at 


Capes Fear, Lookout, Hatteras, and Henry) obtaired from NOS. Comparison 
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water level and wind data at various model lecations were obtained from a CE 
letter report on Hurricane Donna dated 28 April 1961. 

26. Meteorological data for the March 1962 storm were obtained from a 
data base constructed for the WIS project at WES (Resio, Vincent, and Corson 
1982). Comparison water levels and wind data were obtained from a CE letter 


report on the March 1962 storm dated 6 September 1962. 
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PART IV: NUMERICAL MODEL CALIBRATION 
Procedure 


27. The three models in this study were calibrated against the tides of 
the Carolina coast. The tide along the Outer Banks is semidiurnal in nature, 
with the M, constituent predominating. To calibrate the computational models, 
four parameters were adjusted in WIFM: (ay the depths or elevations assigned 
to each cell, (b) the Chezy friction coefficients for each cell, (c) the 
choice of cells used as marigram stations, and (d) the boundary condition 
imposed on the grids. The water depth or land elevation of a cell was esti- 
mated from maps, and only a few deptns were changed slightly during calibra- 
tion. Except in shallow water ana at flow constrictions, the models’ hydre- 
dynamics vere fairly insensitive to changes in Chezy coefficients. The choice 
of which cells to serve as marigram stations was sometimes arbitrary due to 
shoreline approximations and the size of cells relative to hydrodynamic 
details in the prototype. For a model that was acjusted to the point where 
only one WIFM gage disagreed slightly with prototype data, this discrepancy 
was removed by changing the gage's placement in the grid. 

28. Calibrations largely involved the development of proper boundary 
conditions for the models. For @ tide, seawara boundary conditions were 
estimated by wave speed calculations and shoaling factors. With initial 
estimates, the development cf correct boundary conditions proceeded through 
three steps: (a) computations with WIFM using the latest estimate, (b) com- 
parison of computed marigrams to prototype data, and (c) refinement of the 
boundary conditions. Successive iterations vere performed to match computed 
tides against prototype data. The other parameters in WIFM were also adjusted 
during the refinement of the boundary conditions. 

29. The model's ability to accurately reproduce prototype data was 
limited by cell sizes used in the grid (which governed how well topography 
was simulated) and by any approximations used for the boundary conditions. In 
this study, the prototype tides were approximated with just the M, constituent. 
For the North Carolina coast, this approximation provided a good representa- 
tion of the entire tide wnile keeping the expenses of calibration to a mini- 
mum. The number of iterations needed in the boundary condition refinement 


process soars when multiconstituent tides are used, because of the increased 
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number of variables that must be adjusted. 
Offshore Model 


30. The tides of two time periods were simulated in the calibration of 
the offshore model: 5-1U March 1962 (the March 1962 northeaster) and 11 and 
12 September 1960 (Hurricane Donna). Predicted highs and lows for several 
stations along the coast were taken from U. S. Coast and Geodetic Survey 
(1962) tide tables for these two times to serve as prototype data. Four 
stations (Capes Fear, Lookout, Hatteras, and Henry) yielded the initial esti- 
mates for boundary conditions and served as the primary comparisons with WIFM 
computations during the calibration process. Estimates of the M, amplitudes 
and epochs for the grid boundaries were chosea so that the computed one- 
constituent tide would represent the entire predicted tide, and not just the 
My portion of the signal. In order to easily see the agreewent between WIFM 
computations and the prototype data, computed marigrams were ptotted against 
spline-fit curves of the predictea highs and lows. 

31. Boundary conditions were first refined to match the predicted tides 
of the March 1962 storm. The beginning of the WIFM computatior. was chosen to 
be midnight Greenwich Mean Time (GMT) on 5 March 1962. A 3-min time-step was 
used in WIFM, and 111 hr of the tide were computed. Plate i snows computed 
marigrams versus prototype curves for the four Cape stations. There are small 
variances between the peaks of the computed and prototype tides due to the 
lack of other constituents in the grid boundary conditions, but the two tides 
match well for these gages. To check that the model correctly represents the 
prototype tide at other points along the Outer Banks, cther stations are 
compared with computations (Plate 2). Virginia Beach shows an excellent 
match, while the calculated amplitudes at the Currituck Beach Lighthouse are 
about 1/2 ft higher than predicted. This increase is caused by the shoreline 
approximation made in the grid for the barrier island at Currituck, and by the 
proximity of the WIFM gage to this shore. 

32. Plates 2 and 3 illustrate how the crude approximations of the three 
inlets in the offshore grid (each inlet is only one cell wide) induce some 
local amplification of the computed tide. However, the basic hydrodynamic 
Properties of tne inlets are preserved in this model because the cross- 


sectional areas uf the inlets match observations (Jarrett 1976). Tidal wave 
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energy does enter the bays properly, and the marigrams calculated within 
Pamlico and Albemarle Sounds (Plate 3) agree with the information in the tide 
tables. While the NOAA tide tables do not list highs and tows for the bays, 
they note that "In Albemarle and Pamlico Sounds, except near the in_ets, the 
periodic tide has a mean range less than one-half foot" (p. 226, 1962 Tide 
Tables). 

33. After calibration of the offshore model with the March 1962 tide, 
boundary conditions were developed to simulate the predicted tides tor Hur- 
ricane Donna. The starting time of this simulation was noon GMT on il Septem- 
be: 1960, and WIFM computed 33 hr cf tides with a 3-min time-step. The bound- 
ary refinement process was cycled through again, but no minor parareters were 
changed in WIFM. Piate 4 shew: that the computed ampiitudes at Cape Fear fall 
about 1 ft below predictions. This result comes as a compromise: when 
boundary amplitudes were increased in the offshore model near Cape Fear to 
raise the computed peaks, the tid? range throughout the modei became too larg: 
to match predictions at othex stations. Pla-es 5 and 4 compare computafions 
with predictions fer other stations along ths Outer Banks and illustrate the 
tides within the bays. As with the March 1902 tide, these plates show agree- 
ment at Virginia Beach and Curricuck, local increises at the inlets, and the 


proper tides within the bays. 


Nearshore Model 
34. The calibration of the nearshore model was e:ccmplished using 
constituent tide data resulting from harmonic aveiysis o. gage records at 
14 stations within the grid coverage. Although no data collection was per- 
formed for this study, data from ;revious collection efforts were available at 
10 stations in the vicinity of Oregon Irlet, and these were supplemented by 
NOAA data at four additional open-coast locations. Because of its predomi- 
nance, the M, 


2 (lunar-semidiurnal) constituent was selected as a representative 


parameter for model calibration. A summary cf gage locations with M, ampli- 


tudes and epochs is presented in Table 1. Gage locations are Fee in 
Figures 7 and 8. 

35. Initial efforts at calibrating the nearshore model were performed 
using a mean tide from January 1981. First estimates for bathymetry and open 


boundary adjustments were determined using this simulation period. Grid 
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Figure 7. Offshore and nearshore model gage locations 
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resolution in the Oregon Inlet vicinity was modified to relax stability 
exigencies, economize simulations, and yet preserve major inlet features. 
Final adjustments of boundary conditions and friction characteristics were 
performed using a tidal period of 19-22 May 1975 which is coincident with data 
from the previous physical model study (Hollyfield, McCoy, and Seabergh 1983) 
of Oregon Inlet. 

36. Marigrams presenting model versus prototype M, amplitudes appear in 
Plates 7-10 starting at 0300 (GMT) 20 May 1975. The model quite accurately 
reproduces the open-coast tide as observed at all available stations. Model 
performance at Oregon Inlet is characterized by good agreement at Bonner 
Bridge and the Pea Island Coast Guard Station. The remaining stations present 
comparisons at bayside locations in the inlet vicinity. They also show good 
agreement with amplitude variances iess than Q.1 ft and some slight phase 


shifting attributed to channel and shoal features subscale to grid resolution. 
Shore Process Model 


37. Unlike the coarser models used in the study, the shore process 
model presented a unique case tor calibcation. Open-boundary conditions were 
not adjustable since they were transferred from simulations on the nearshore 
model. Because the primary function of this model was to provide hydrodynamics 
for sedimer.t transport studies, cell depth changes were restricted in order 
to maintain the same bathymetry for these two modeling efforts. Friction 
values were adjusted to coincide with representations in the nearshore model. 
The inclusion of the momentum advection terms and the very fine grid resolution 


(Ax = 100 ft) with the resultant small time-step (At < 15 sec) required 


ee, adjustments of the eddy viscosity coefficient to achieve simulation 
stability, and yet produce realistic circuiation currents and horizontal 
eddies evident at the inlet under normal and proposed configurations. Cali- 
bration efforts for this medel were performed for mean tides produced by 
nearshore model simulations in the aforementioned May 1975 period. Results 


are discussed in PAR V, and model gage locations are shown in Figure 38. 
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PART V: NUMERICAL MODEL VERIFICATION 


38. After calibration, the ability of the models to simulate the 
nydrodynamics of storm surges was verified with the two historical storms. 
Hurricane Donna and the March 1962 northeaster produced some ot the highest 
surges recorded along the Outer Banks, and they had severe effects on Oregon 
Inlet. These storms were considered rigorous tests for the models of the 
North Carolina coast. 

39. The eye of Hurricane Donna passed just west of Pamlico Sound 
during the storm's movement north (Figure 9), with the highest winds concen- 
trated over the Sound. During the peak of the storm, 60- to 80-knot winds 
blew across Pamiice Sound from the south and southwest, causing a large surge 
to pile up along the bay side of the barrier islanus near Oregon Inlet. The 
elevation difference between the surge in the bay and the level of the Atlantic 
Ocean created a 10-ft head across the inlet, which caused an enormous flow 
through the inlet. The surge in the bay also caused extensive flooding at the 
barrier islands, Roanoke Island, and the nearby mainland. 

20. The March 1962 storm derived its force not from high peak winds but 
from sustained gales over several days. The iow pressure system which drove 
the storm remained along the eastern seaboard for over a week, generating 35- 
to 50-knot winds from the north and northeast from 6 March to 8 March. These 
gales caused a surge on the Atlantic side of the Outer Banks near Oregon 
Inlet and a drawdown on the east side of Pamlico Sound. The head difference 
at the inlet caused a large flow into the bay through an entire tidal cycle on 
7 March, which, together with setup, created flooding in nearby Roanoke Island 
and the mainland. 

41. The storm simulations used the historicai tides developed during 
calibration as initial and boundary conditions for the nearshore and offshore 
models. The initial conditions for a WIFH storm run were supplied by data 
sets (dubbed "hotstarts") saved during tidal computations. Hotstarts contain 
the free-surface elevations and horizontal velocities at every point in a 
grid, which enables one to start W!IFM calculations with this information as a 
realistic, dynamic initial condition. The same tidal boundary conditions were 
also used during storm simulations, so that the tide appears throughout each 
storm hydrograph. Hence, computed hydrographs can be directly compared with 


historical marigrams. 
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Figure 9. Track of Hurricane Donna 


42. With the boundary and initial conditions for a storm simulation 
specified by tides, the meteorological forces appear in the computations as 
the terms ES 5 Lee » and 5 in Equations 2 and 3. These terms represent the 
shear (wind stress) and normal (barometric pressure differential) forces 
applied to the water surface during the storm. WIFM obtains values for the 
wind and pressure fields from datasets created by other computer codes, such 
as the SPH Program. Wind velocity and pressure head values for every grid 
cell are read into WIFM every few time-steps, and then WIFM converts the wind 
velocity into surface shear stress by Charnock's linear relation (Garrett 


1977). The drag coefficient for Charnock's method is: 


to 
co 
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Cy = 0.00075 + 0.000067W (19) 
where W is the 10-metre wind speed in metres per second. The surface 
stress is calculated by Taylor's (1916) equation: 

T= PPG w (20) 
D 

where p is the air density, assumed constant for this study. The total 
surface stress then is resolved into components for the ae and a terms in 
WIFM. 

43. Due to the increase in flow velocities through Oregon Inlet during 
the storms, the computational time-step in WIFM had to be adjusted for the 
nearshore model. Both Hurricane Donna and the March 1962 northeaster were 
simulated with a 3-min time-step in the offshore model (the same as for the 
tide runs). The March 1962 computation begins at midnight GMT on 5 March and 
goes to noon GMT on 9 March. The Hurricane Donna simulation begins at noon 
GMT on 11 September 1960 and computes to 9 p.m. GMT on 12 September. For the 
nearshore model, the time-step was reduced to 50 sec but the same starting and 
ending times were kept. The March 19€2 storm was simuiated with a 60-sec 
time-step from hours 33 to 78 (9 a.m. GMT on 6 March to 6 a.m. GMT on 8 March). 
The simulation of only 45 hr of this storm greatly reduced computational 
costs. 

44. Two letter reports of SAW (Davidson 1961, Grygiel 1962) contain 
the prototype storm surge data used in this study. Surge histories for the 
peak pericds (7 March 1962 and 11 and 12 September 1960) for most of the sites 
in Figure 7 can be found in these reports. No prototype data were available 
at or near Oregon Inlet, but there are recorded marigrams available from 
Rodanthe, Nags Head, Stumpy Point, and Point Harbor. These four sites all lie 
within 25 miles of the inlet, so that agreement between them and calculations 
can indicate the accuracy of calculations at Oregon Inlet. Tabulated measure- 
ments of the tide highs and lows for the months of March 1962 and September 
1960 provided independent checks on the accuracy of the surge histories; these 
two sources agreed completely. The tabulated data also gave an indication as 
to where a "mean water level" may lie for each prototype gage. An average of 
the month's highs and lows (excluding the storm period) seldom came close to 


the datum listed with the tabulations. Differences are mainly attributable 


to two factors: datum discrepancies and local superelevation. The datum for 
all tide gages was based on the National Geodetic Vertical Datum (NGVD), which 
is not intended to equal the local mean water level. Bench marks used tc 
reference the gaces were not tied together by a single survey and hence ..ay 
not lie in the same datum plane. In addition, local mean water levels, par- 
ticularly in the sounds, should be above NGVD due to a general rise in sea 
level, since NGVD was established in 1929, and also due to local supereleva- 
tion of the bays caused by freshwater inflow into the sounds. Treshwater 
inflow, and its local effects, was not considered in this modeling effort. In 
order to account for the above factors, datum adjustments were estimated as 
differences between the NGVD values and the monthly high/low averages. Table 3 
summarizes the sites used in the verification of the two largest models, gages 
within the grids, and the datum adjustments estimated from the high/low 
averages. Tabulated data were not available for all of the prototype mari- 
grams, and it must be noted that the March 1962 surge at Nags Head is esti-~ 
mated since the tide gage there was lost during the storm. Grygiel states 
that the estimate "is based on reliable evidence and should be less than 1 ft 


in error..." (p. 3). 
Offshore Model 


45. Plates 11-13 compare the computed storm surge cf the March 1962 
storm with the prototype data for the 11 sites listed in Table 3. Where datum 
adjustments are available, they have been subtracted from the prototype 
records in these plates. While the comparisons agree well, some differences 
occur due to the shoreline approximations used in the offshore model and due 
to the presence of wave setup in the prototype data. Computations for 
OrientatL, Rodanthe, Point Harbor, and Columbia exhibit excellent fits to the 
field records. Minnesott Beach and Cherry Point have computed peak suiges 
that are slightly high due the shoreline approximations made for the Neuse 
River. The cell sizes in the offshore grid necr Stumpy Point Bay also do not 
allow this small bay to be well represented; therefore the computations and 
prototype for this site differ. The presence of wind-induced wave setup in 
the prototype records of Wrightsville Beach, Hatteras, and Englekard can 
explain why these computations fall below field measurements. The divergence 


in elevations after hour 60 (the peak of the storm) at Hatteras shows how wave 
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setup can act to keep the water elevation raised in the prototype data. Since 
WIFM cannot model this effect, the computed water level at Hatteras begins 

to fall after the peak of the stcim. Prototype data for Nags Head are not 
available (the tide gage disappeared during the storm) but estimates of 

the surge have been made by SAW from other information. shis estimate ap- 
pears with computations for Nags Head in Plate 13 with no datum shift esti- 
Mates removed. The computed and estimated prototype records agree in form, 
and the roughly 2-ft elevation difference can be accounted for by datum 
shifts, wave setup, and errors in Grygiel's estimate of the record at Nags 
Head. 

46. Plates 14-16 compare the computed storm surge of Hurricane Donna 
with prototype data for the offshore model. Ten sites are shown, and the 
datum shift estimates of Tabie 3 have been removed from the prototype 
records. Atlantic Beach and Nags Head, both sites on the open ocean, show 
the best agreement. For sites in the bay, the effects of wave setup and 
shoreline approximations play an important role in the agreement between 
computations and field records. An added factor is the speed at which Hur- 
ricane Donna moved across the North Carolina ceast. The eye of the storm 
traversed the distance from Cape Tear to Cape Henry in about 8 hr, causing 
rapid changes in wind speed and direction for sites in Pamlico and Albemarle 
Sounds. The rapid movement of the storm, along with the shoreline simpli- 
fications of the Neuse and Pasquotank Rivers, caused some phase difference 
between computations ard the prototype records at Cherry Point, Minnesott 
Beach, and Elizabeth City. Columbia, a site up a narrow recess of Albemarle 
Sound, also has a phase shift. The prototype records for Hatteras and 
Nags Head show the effect of wave setup at the peak of Hurricane Donna 
(hour 22) but these records agree well with computations. Despite the shore- 
line approximations for Stumpy Point Bay in the offshore grid, the results 
match well at this site. The computations at Englehard overpredict the 
surs,e, but the results agree at hour 22, the time of maximum flow at Oregon 
Inlet. 

47. Plates 11-16 illustrate that the offshore model can simulate both 
historical northeasters and hurricanes to a reasonable accuracy, with a mini- 
mal number of grid cells. Considering the accuracy with which the datum 
shift adjustments in Table 3 are estimated, the comparisons of computations 


with field data for the offshore model are quite good. 
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Nearshore Model 


48. Plates 17-19 illustrate the excellent agreement between compu- 
tations for the March 1962 storm in the nearshore model and the prototype 
data. Once again, the datum shifts listed in Table 3 are included in the 
field records. The computations match all 11 sites quite closely. Davis, a 
gage not represented in the offshore model due to cell sizes, kas no datum 
shift included in the prototype record. It can be expected that any shift 
would be of the same order as nearby sites such as Wrightsville Beach. The 
computations of Stumpy Point agree with the prototype because Stumpy Point Bay 
can be reasonably approximated by the small cells in the nearshore grid. 

49. Plate 20 shows the flow patterns through Oregon Inlet for the peak 
of the March 1962 northeaster (hour 60, or noon GMT on 7 March) in the near- 
shore model. Piates 21 and 22 show marigrams for eight locations (Figure 3) 
in Oregon Inlet, while Plates 23-26 show the corresponding velocity records. 
The velocities at the inlet tell what happened during the storm. The normal 
ebb-flood cycle of tidal flow through the inlet was halted by the surge sea- 
ward of the inlet on 7 March. For a 22-hr period »etween 8 a.m. GMT on 
7 March (hour 56 in the calculations) and 6:00 a.m. GMT on 8 March (hour 78), 
the flow was entirely into Pamlico Sound. According to computations, the peak 
flow velocity was about 9 fps, or about twice the normal tidal velocity. At 
the peak of the storm about 630,000 cfs of water flowed into Pamlico Sound. 
Figure 10 shows the computed surge contours at noon GMT on 7 March (hour 
60) for the area of Oregon Inlet, Roanoke Island, and the nearby mainland. 

The head difference driving the flow is about 5 to 6 ft. 

50. Comparison of the nearshore model calculations of Hurricane Donna 
with the prototype data can be found in Plates 27-29. Most of the phase 
differences found in the offshore model computations do not appear “n the 
nearshore model due to the increased resolution in this grid. As with the 
offshore model, the effects of wave setup appear in the plots for Hatteras and 
Nags Head. Considering the severity of Hurricane Donna, and the resolution of 
the grid for the nearshore model, Plates 27-29 show that this model can 
simulate a historical hurricane. 

51. Plate 30 illustrates the severity of the ebb flow through Oregon 
Inlet at the peak of Hurricane Donna (10:00 a.m. GMT, 12 September, or hour 22 


of the computations). Plates 31 and 32 show the surge histories of the 
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locations in Figure 9, while Plates 33-36 show the corresponding velocities. 
Once again, the velocity records describe how the storm affected the inlet. 
The buildup of surge in Pamlico Sound on the bayward side of the barrier 
island caused by the high winds from the southwest drove a violent outflow 
through the inlet between 7:00 a.m. and 4:00 p.m. GMT on 12 September. Fig- 
ure ll shows the computed contours for the area around Oregon Inlet at the 
peak of Hurricane Donna, and this figure illustrates the 10-ft head difference 
across the inlet which caused the ebb flow. Peak velocities of 12 to 14 fps 
occurred in the inlet, with 18 fps calculated at the north oceanside approach. 
The nearshore model indicates that about 1 million cfs of water flowed out 


through Oregon Inlet at the peak of Hurricane Donna. 
Shore Process Model 


52. Prototype storm data were unavailable for the relatively small 
region encompassed by the shore process model. Since it was employed princi- 
pally for tidal simulations, results were compared for the My constituent 
tide. Plate 37 shows marigrams fur a 16-hr period starting at 1900 (GMT) 

20 May 1975. The model produced a good match of elevations with prototype 
data at locations just seaward, within, and west uf the inlet, with variances 
of tess than 0.1 ft, and with slight phase shifts on the bayside due to 
features subscale to this grid, and to bayside -oundary conditions affectea by 
similar scale restrictions in the uriving nearshore model. Velocity focords 
(Plate 38) in the inlet and Davis Slough were reproduced with proper phasing 
of the peak flood and evb flows and were generally within 0.5 fps of observed 
data indicating proper tidal exchange. Tidal circulation patterns matched 
quite well in direction and magnitudes with reported data of the same period 


as presented by Hollyfield, McCoy, and Seabergh (1983). 
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PART VI: CONTROL STRUCTURES IMPACT 


53. The Oregon Inlet control structures would consist of two parallel 
jetties extending seaward from the ends of Bodie and Pea Islands. Jetty 
spacings of 2,500, 3,500, and 5,000 ft (Figure 12) were evaluated in both the 
nearshore and shcre processes model. In the nearshore model, both normal 
astronomical tides and an extreme meteorological event represented by Hurricane 
Donna (September 1960) were simulated. The shore processes model was used to 
develop much finer flow detail around Oregon Inlet for eventual input to the 
sediment transport simulations. In order to dztermine the maximum influence 
changes at Oregon Inlet could have on bay-side storm tides, tests were run 
with Oregon Inlet closed. These numerical simulations provided good estimates 
of how the proposed structures at Oregon Inlet would affect surrounding areas 


during both normal and storm generated tides. 
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Figure 12. Proposed alternate jetty configurations 
for Oregon Inlet 


Preliminary Tests 


54. The eastern edges of the offshore and nearshore grids (Figure 2) 
were originally placed far enough from Oregon Inlet so that any changes to the 
inlet would not disturb the boundary values of the models. This assumption 
was checked with the offshore model by closing the iniet (the severest change) 
and then simulating Hurricane Donna (the severest surges). Results of this 


run were then compared with a similar run with the inlet open. The marigrams 
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at the nearshore grid boundary east of Oregon inlet agreed completely, indi- 
cating that this nearshore boundary was indeed far enough from the inlet. 

55. Hurricane Donna also was simulated in the nearshore grid with the 
inlet closed. This test served to answer the question "What is the maximum 
radius of influence at which Oregon Inlet can affect a peak storm surge?" A 
comparison of nearshore computations with existing conditions with the closed 
inlet case illustrates the effects of closure. Plates 39 and 40 show the two 
cases for the sites in Figure 8. At the closed inlet the surge differences 
are, of course, large; but they diminish rapidly away from the opening. The 
difference is only 1 ft at the Roanoke Sound channel gage, a site 3 miles from 
the center of the inlet. Tigure 13 shows the peak surge contours around 
Oregon Inlet for the closed inlet and for natural conditions. The surge 
differences become nil near Rodanthe south of the inlet (13 miles), Manteo 
north of inlet (12 miles), and 10 miles west of the inlet. So, the closed- 
inlet computation with the nearshore model indicates that the maximum radius 
of influence cf the inlet is about 10 miles but changes of a foot or more only 
will occur within a radius of 4 miles from the inlet's center. For example, 
water elevations at the Pea Island Coast Guard Station and the Oregon Inlet 
Marina increased by 2.4 ft and 1.4 ft, respectively, under closed inlet con- 
ditions. Any ronradical changes to the inlet's hydrodynamics (such as 
jetties) will exert a minimal influence on peak surge levels, as the jetty 


computations of this study will show. 


Storm Surge Tests 


56. All of the storm surge simulations with the jetties were computed 
in the nearshore model. Hurricane Donna was simulated with the 2,500- and 
5,000-ft structures. The 5,000-ft structure changed the surge levels and peak 
velocities very little from the nearshore simulation with existing conditions 
at the inlet, while the 2,500-ft structure slightly raised the surges and 
velocities at the peak of Hurricane Donna. 

57. Plate 41 shows the flcw patterns through Oregon Inlet at the peak 
(hour 22) of the high velocities that occurs nearest the north jetty, due to 
the presence of shallow water near the south structure. Plates 42 and 43 com- 
pare marigrams of the 2,500-ft structure to existing conditions for the sites 


of Figure 8 while Plates 44-47 compare the corresponding velocity records. The 
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largest change in the peak elevations occurs at the North Ocean Entrance, a 
site at the ends of the jetties. The constricting effects of the structure 
raise the peak surge avout 2 ft in the middle cf the inlet, while little 
difference is noted on the bay side of the opening. Velocities on the bay 
side are actualiv reduced by the jetties (i.e. Bonner Bridge) until the flood 
flow reaches the middle of the inlet. Figure 14 shows a comparison of eleva- 
tion contours around the inlet for the peak of the storm for this structural 
alternative. 

58. Piate 48 shows the flow patterns at the inlet for the 5,600-ft 
structure. The shoal in the channel acts to move the rastest velocities near 
the two jetties. Plates 49 and 50 compare water elevations for this case with 
existing conditions, while Plates 51-54 compare velocities for the stcrm. The 
influence of the jetties on either the surge or the velocities is negligible, 
even in the channel. Figure 15 shows a comparison of computed surges at the 
peak of the storm, and clearly indicates that the effects of the 5,000-ft 
spacing jetties are limited to small variances within the inlet, and no 
impact on Pamlico Sound, as represented by the nearshore model. 

59. A tabular comparison of computed peak surge results for Hurricane 
Donna at selected Outer Bank locations with the various inlet configurations 
is presented in Table 4. Results in this comparison also indicate no impact 
of the structures on peak surge levels within the Albemare/Pamlico Sound 


systems as recorded along the Outer Banks. 


Tidal Tests 


60. Tidal simulations with the jetties were conducted with the near- 
shore and shore process models. A mean tide was simulated with the 2,500-, 
3,500-, and 5,000-ft jetties. Structure effects at the inlet vicinity were 
evaluated with the shore process model, and effects at more distant locations 
along the Outer Banks were determined with the nearshore model. The jetty 
configurations were simulated with natura] bathymetry features including the 
offshore shoal and no specific dredging. In general, structure effects were 
limited to the immediate inlet vicinity with no changes observed at Outer Bank 
locations. 

61. Plates 55-59 present comparisons of natural conditions versus the 


2,500-ft jetty alternative for both nearshore and shore process models. Local 
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Figure 14. Comparison of computed surg2 contours at peak of 


Hurricane Donna (natural inlet versus 2,500-ft jetty spacing) 
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Figure 15. Comparison of computed surge contours at peak of 
Hurricane Donna (natural inlet versus 5,000-ft jetty spacing) 


4). 


variations of velocity and elevation were apparent in the inlet and within the 
jetties. The Pea Island Coast Guard Station experienced range reductions of 
0.30 ft while the Davis Slcugh location experienced a lesser reduction of 

0.10 ft. At stations along the Outer Banks such as Bodie Island Coast Guard 
Station, Nags Head or Rodanthe, amplitude variations were negligible or non- 
existent. Stations located farther within the Sound showed no variation. 

62. Similar comparisons of natural conditions versus the 3,500-ft jetty 
are shown in Plates 60-64, while Plates 65-69 show the 5,000-ft jetty results. 
As anticipated, these larger jetty spacings produced smaller variances at the 
inlet proper, and no discernible changes in range along the outer banks or 
within the Sound. Plates 70-75 show jetty alignments and flow patterns in the 


nearshore model at periods of peak flow and ebb tidal flows. 
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PART VII: SHORE PROCESS MCDEL REQUIREMENTS 


63. Concurrent with the numerical model efforts described in this re- 
port was a study of sediment transport under current and wave-induced current 
interactions at Oregon Inlet. The application of the WIFM code to the common 
computational grid used by both studies has been referred to as the shore 
process model in this report. WIFM was employed to provide elevation and 
velocity data at ])-min intervals for each celi in the computational grid to 
the sediment transport models. Suck data were required for the following; 

a mean, spring, and neap tide, a mean tide with 2,500, 3,500-, and 5,000-ft 
jetty spacings, and south jetty only; and for the March 1962 northeaster. 

64. An M, constituent tide was again chosen as a boundary condition. 
The My amplitude was adjusted to represent the full tide values for mean, 
spring, and neap conditions as reported at Nags Head, North Carolina, by NOAA 
Tide Tables. Am»olitudes of 1.6 ft, 1.9 ft, and 1.3 ft were chosen for the 
mean, spring, and neap tides, respectively. Appropriate boundary conditions 
were developed for the nearshore medel to produce the proper tides and to 
provide connecting boundary conditions ror the shore process model. Marigrams 
at typical statioas from shore process model simulaticns under mean tide 
conditions are presented in Plates 76 and 77. Flow patterns at flood, slack 
after flood, and ebb portions of the mean tidal cycle at the inlet are shown 
in Plates 78-80. They indicate the model's ability to realistically produce 
the rather complicated horizcntal flow patterns through the inlet as affected 
by channelization and shoaling. Similar results for spring tides are pre- 
sented in Plates 81-85 while neap tide conditions appear in Piates 86-90. 

65. The sediment transport studies alse required hydrodynamic informa- 
tion for the various structural configurations at the inlet. These alterna- 
tives were simulated using the mean tide condition. Simulations with the 
2,500-, 3,500-, and 5,000-ft jetty spacings with the nearshore model were used 
to provide driving boundaries for equivalent cases on the shore process model. 
Marigrams showing typical shore process model results for these cases are 
presented in the plates describing structural effects. Flow patterns with the 
jetties in rlace for various portions of the tidal cycle appear in Plates 91l- 
99. Flow redistribution around and through the jetties is clearly depicted as 
well as the Jecal acceleration near the offshore shoal. 


66. A fourth structural alternative was simulated for the sediment 
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transport studies. This was a single structure consisting of the south jetty 
only from the 5,000-ft spacing alignment. A comparison of the 5,000-ft jetty 
effects and existing conditions (Plates 100 and 101) demonstrated that the 
driving boundary conditions from the nearshore model were unaffected by the 
jetties; therefore 2 nonstructural mean tide boundary condition from the 
nearshore model was applied to the shore process model. Marigrams from this 
simulation appear in Plates 102-104 with the fiow patterns presented in 
Plates 105-107. 

67. For sediment transport studies under severe conditions, 24 hr of 
the March 1962 northeaster were simulated starting at 0600 (GMT) on 7 March. 
Marigrams at typical stations appear in Plates 108-110 with flow patterns near 
the peak of the stcrm at noon shown in Plate 111. The extreme flood flows at 
the inlet are auite apparent. As indicated in Plate 110, flc* through che 
inlet was unidirectional into the bay for the entire simulation period as 


observed during the actual ctorm. 
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PART VIII: SUMMARY AND CONCLUSIONS 


68. Numerical hydrodynamic models were developed ana tested for the 
purposes of evaJuating the influence of proposed structures under storm con- 
ditions and providing elevation and velocity data for concurrent numerical 
sediment transport studies at Oregon Inlet, North Carolina. Data collected 
for previous physical model studies were supplemented by additional data from 
NOAA and unpublished SAW letter reports for use in calibrating and verifying 
tne models under existing tidal conditions and for two severe storms of 
record--the March 1962 northeaster and Hurricane Donna (1960). Three nodels 
with progressively finer resolution (offshore. nearshore, and shore process) 
were successfully calibrated (relative to their respective degree of resolu- 
tion) to replicate an astronomical tidai event. Using a tuned hurricane 
windfield model for Donna and windfields developed by the Wave Information 
Study at WES for the March 1962 storm as forcing parameters, the models were 
able to duplicate observed marigrams throughout the model area. Observed bay- 
ocean head differentials also were modeled correctly, particularly the 20-hr 
plus flood event through Oregon Inlet during the peak of the March 1962 storm. 

69, Two structural alternatives (involving parallel jetties with 
2,500- and 5,000-ft-wide spacings) were studied under tidal and storm condi- 
tions. Structural effects were determined to be limited to the inlet under 
tidal conditions and to the immediate inlet vicinity under cixcumstances 
approximating Hurricane Doura. The maximum possible influence exerted by anv 
changes at the inlet was determined by simulaticns with the inlet completely 
closed (this was not actually a proposed improvement plen) under Hurricane 
Donna conditions. For this case, no change was noted beyond a 12-mile radial 
distance from the inlet. During Hurricane Donna, peak surge levels at the Pea 
Island Coast Guard Station and the Oregon Inlet Marina were increased by 2.4 ft 
and 1.4 ft, respectively, with total closure cf Oregon Inlet. These stations 
are in the immediate inlet vicinity (within a 1.2-mile radius). For the 
2,500- and 5,000-ft spacings, peak surge level increases for the nearby sta- 
tions mentioned above were even smaller, with little or no differences noted 
at farther distances from the inlet. 

70. Finally, simulations at very fine grid resolutions were made to 
provide hydiodynamic data for the sediment transport studies covering the 
normal range of tides with no jetties, for four structural alternatives with 
a mean tide, and for the historical March 1962 northeascter with no jetties. 
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Table 4 


Comparison of Hurricane Donna Peak Suse Levels Along Outer Banks 


Distance from 2,500-ft 5,000-ft 
Oregon Inlet, Jetty Jetty 
_ Location (Sound Side) miles _—s—s—« Existing Closed Spacing Spacing 
Bodie Island: 
Oregon Inlet Marina abe 7.8 9.2 8.1 7.8 
Bodie Island 3.0 Thee) 8.8 8.1 7.9 
Lighthouse 
Headquarters Island 8.5 6.1 6.4 6.2 6.1 
Sound Side (Nags Head) 13.25 6.6 6.8 6.7 6.6 
Pea-Hatteras Island: 
U. S. Coast Guard i.0 Coy) 38.9 Hal! 6.5 
Station 
Goose Island 4.5 8.4 8.9 8.5 8.4 
U. S. Fish & Wildlife 7.6 8:2 8.4 8.2 8.1 
Service Headquarters 
Rodanthe 13.6 7.8 7.8 As8 7.8 
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APPENDIX A: NOTATION 


Regional constants derived from stretching transformation of 
coordinate system 


Coefficient matrices 

Chezy frictien coefficient 

Wind drag coefficient 

Tctal depth of water column d = 7-h 

Coriolis parameter 

Terms representing external forces such as wind stress 
Acceleration due to gravity 

Local ground (cell) elevation above datum 

Integer time-step counter 


Dimensionaless parameter used to characterize stability 
criterion 


Rate of water volume change in the system (for example 
through rainfall or evaporation) 


Time 

Vertically averaged water velocity in x-direction 
Matrix consisting of mn, u, and v as functions of x, y andt 
Vertically averaged water velocity in y-direction 

Largest velocity encountered at a computational cell 

10-metre wind speed 

Cartesian coordinate system axes names 

Smaller value of Ax and dy 

Time-step 

Length of computational cell in x- and y-directions 

Eddy viscosity coefficient 

Water-surface elevation above datum 


Hydrostatic water elevation due to atmospheric pressure 
differences 


Two-dimensional differences operators 
Air density 
Surface stress of wind 


Intermediate time-step level 


Positive integer representing computational grid line 
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